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Searching for CU Vir-type cyclotron maser from a Ori E: 
The role of the magnetic quadrupole component 
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ABSTRACT 

In this paper we present new and archive radio measurements obtained with the 
Very Large Array of the magnetic chemicahy peculiar (MCP) star crOriE. The radio 
data have been obtained at different frequencies and are well distributed along the 
rotational phases. We analyze in detail the radio emission from a Ori E with the aim to 
search evidence of circularly polarized radio pulses. Up to now, among the MCP stars 
only CUVirginis shows 100% polarized time-stable radio pulses, explained as highly 
directive electron cyclotron maser emission, visible from Earth at particular rotational 
phases, like a pulsar. Our analysis shows that there is no hint of coherent emission at 
frequencies below 15 GHz. We conclude that the presence of a quadrupolar component 
of the magnetic field, dominant within few stellar radii from the star, where the maser 
emission should be generated, inhibits the onset of the cyclotron maser instability in 
crOriE. 

Key words: masers - stars: chemically peculiar - stars: individual: aOriE - stars: 
magnetic field - radio continuum: stars. 



1 INTRODUCTION 

Magnetic chemically peculiar (MCP) stars present photo- 
metric, spectroscopic and magnetic variability with a sin- 
gle period. As to the magnetic field, the most impor- 
tant observable is the efi^ective or longitudinal magnetic 
field (-Be), that is the average over the stellar visible 
disk of the longitudinal components of the magnetic field 
IjLeone. Catalano &: Ca tanzaro 2000'). To ex plain t h e pre - 
vious phenomenology [Babcock (1949|) and IStibbd (llQSd ) 
proposed the presence of a magnetic dipole, whose axis 
is tilted with respect to the rotational axis, the so called 
Oblique Rotator Model, and a non-uniform distribution of 
some chemical elements on the stellar surface which ro- 
tates rigidly. Spectral, light and magnetic variability would 
be a consequence of stellar rotation. Later, it has been 
cleared tha t the MCP stars can present multi polar mag- 
netic fields. iBvchkov. Bvchkova fc Madeil (120051 ) collect all 
available MCP magnetic curves and they show that in many 
cases their shape is more complex then a simple sinusoidal 
wave. 

Non-thermal ra dio emission is observed fr om about 25% 
of the MCP stars JLeone. Trigilio fc Umana| [i994'). In ac- 
cord with the oblique rotator model, the radio emis sion is 
also v ariable as a consequence of the stellar rotation (JLeonel 
[l99l|), suggesting that the radio emission arises from a stable 
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optically thick co-rotating magnetosphere. Radio emission is 
ascribed to a radiatively-driven stellar wind. The gas flow 
brakes the magnetic field lines far from the star (Alfven sur- 
face) forming current sheets, where electrons are accelerated 
to the mildly relativistic regime. Energetic electrons propa- 
gate back, along thin magnetospheric layers named middle 
magnetosphere, to the inner magnetospheric regions radi- 
ating by gyr osynchrotron emission mechanism at the radio 
wavelengths (|Trigilio et al.ll2004 [beto et al.ll2006l 'l. 

Out of MCP stars, CUVirginis (HD 124224 = HR5313) 
is the only known source characterized by broadband, 
highly polarized and time - stable pulses at 1 . 4 and 2.5 GHz 
(i Trigilio et aLll2000l.[20og I2OIII : iRavi et al.l I2OI0I : ILo et al.l 
2OI4). IStevens fc Georgg l|2010l ) have also reported a pulse 
detection at 610 MHz. Radio pulses stand out for 1 dex over 
the continuous emission, the single pulse duration ranges 
from 5 to 10% of the rotational period and are observed 
in coincidence with the null values of Be, that is when the 
magnetic dipole axis is perpendicular to the line of sight. 

Such a behavior has been ascribed to the electron cy- 
clotron maser (ECM) emission mechanism powered by an 
anisotropic pitch-angle (angle between electron velocity and 
local magnetic field) distribution that the non-thermal elec- 
trons, responsible of the gyrosynchrotron stellar radio emis- 
sion, can develop propagating in a magnetic flux tube to- 
wards regions of increasing magnetic field strength. Elec- 
trons with an initial large pitch-angle are soon reflected out- 
ward due to the magnetic mirroring, whereas electrons with 
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a small pitch angle can reach the inner magnetospheric lay- 
ers where they are thermalized in the dense plasma. The 
reflected electron population is characterized by a pitch an- 
gle distribution deprived by the electrons with a small pitch 
angle. This mechanism amplifies the extraordinary magneto- 
ionic mode producing nearly 100% circularly polarized radi- 
ation at frequencies very close to the first or second harmonic 
of the local gyro-frequency {v-b = 2.8 x 10~^i3/G GHz), in 
a dire ction almost perpendi cular to the local magnetic field 
lines (JMelrose fc DuM 19821 ). However the fundamental har- 
monic is probably suppressed by the gyromagnetic absorp- 
tion. 

The ECM mechanism has been considered to ac- 
count for the strongly polarized, intense an d narrow band 
s hort time spikes ob served in the Su n (IWillson 



Winglee & Dulk"l986'), dMe flare stars ("Lan g et al 
Lang & Willson 1988; Abada-Simon ct al. 199' 
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RSCVn binaries (jSlee. Havnes fc Wrightlll984i r iOsten et all 
120041 ). ECM explains the planetary low frequency emis- 
sion, in particular the Jupiter decametric radiati on and the 
Earth 's Auroral Kilometric Radiation (AKR) (|Treumannl 
l2006l ). 

The fully polarized pulses observed on CUVir are 
broadband and persistent over long timescale (years). The 
large bandwidth is a consequence of the wide range of the 
magnetic field strength in the region where the maser am- 
plification occurs. The observed maser emission is the su- 
perimposition of narrowband emission from different rings 
above the magnetic poles. The continuos supply of non- 
thermal electrons, developing a loss-cone anisotropy, main- 
tains stable the electron cyclotron maser emission. Follow- 
ing th e tangent plane beaming model, proposed for the 
AKR (JMuteL Christopher fc Picketdliooj ) and successfully 
applied to explain the narrow peaks observed on CU Vir 
l|Trigilio et al.ll201H ). the amplified radiation is beamed tan- 
gentially to the polar ring where the cyclotron maser insta- 
bility takes place. Then, during the propagation through the 
denser magnetized plasma of the inner magn etosphere, the 
radiation is refracted upward by a few degrees (JTrigilio et al.l 
|2011I : ILo et al.1120121 ). Since the magnitude of this angle de- 
pends on frequency, the radiation is detected in different mo- 
ments during the rotation of the star, causing a frequency 
drift of the observed pulses. 

Highly polarized broad-band radio pulses, still ex- 
plained as ECM, have been also rec ognized in t he ul- 
tra cool main seque n ce dwarf stars (IBerger et al.l 



tra cool mam seque n ce dwari stars ( | £ierg : er et al 
Burgasser fc PutmanI l2005l : iHaUinan et al.l l2006l . 
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20081 ). Despite the great difference of the physical charac- 
teristics between these two extreme classes of stars, late M 
and MCP stars show similar behavior at the radio wave- 
lengths. The existence of a well ordered and stable axisym- 
metric magnetic fie ld in the rap idly rotating fully convec- 
tiveM type stars (JDonati et al. 2006a, 2008; iMorin et aD 
l2008al la) can be the reason of the observed similarity. 

The study of the radio pulses from CU Vir also provided 
the possibility to evidence variations of t he stellar rota.tional 
perio d with a high degree of confidence (|Trigilio et al.ll2008l . 
uOllh . Discovering the CU Virginis type coherent emission in 
other MCP stars would provide an useful tool for the study 
of the angular momentum evolution in this kind of stars. 

On the basis of the results gathered from CUVir, the 
possibility to observe the same type of coherent emission 



from other MCP stars requires an appropriate stellar mag- 
netic field geometry, with B^ presenting at least a null 
value during the stellar rotation. Among the already known 
MCP stars presenting radio emission o Ori E (HD 37479 
= HR 1932) is an ideal candidate to search for the pres- 
ence of this type of coherent emission. Its magnetic field 
presents the appropriate geometry, being the rotation axis 
inclinatio n i = 72° an d magnetic axis obliquity /3 = 56° 
(Bohlend er et all |l983) . 

The coher ent pulses have been observed from CU Vir 
(Bp = 3000 G. iTrigiho et al.1 l|2000l )) at low frequencies (< 
2.5 GHz). We expect that ECM emission from crOriE could 
be observed at higher frequencies, as a consequence of t he 
stronger magnetic field (Bp = 6800 G. lTrigilio erahl l|2004h ). 
In this paper we present new multifrequency (1.4, 5, 8.4, 15, 
22 and 43 GHz) VLA observations of cr Ori E. In addition, 
to obtain as complete as possible radio light-curves, we have 
retrieved all unpublished VLA archive data. The primary 
aim is to search for amplified emission strongly polarized, 
but we can also probe different layers of the o Ori E mag- 
netosphere, because of the dependence of the frequency of 
the gyrosynchrotron radio emission from the magnetic field 
strength. 



2 a ORI E 

a Ori E is a magnetic helium-strong star of spectral type 
B2V, mass M. = 8. 9 Mq JHu ngcr, Hcbcr & Grootc 198<|), 



radius 7?* = 4.2 R© (IShore fc Brownlll990l) . located at a dis- 
tance of about 350 pc (ESA|[1997|). It is significantly hotter 
then CUVir (spectral type AOV), this involves a stronger 
radiatively driven wind and, therefore, a greater mass loss 
rate. On the other hand, the polar magnetic field strength 
of crOriE is 6800 G, versus 3000 G of CUVir, and the ro- 
tational period is 1.19 , versus 0.52"*. All these two factors 
infiuence the confinement of the wind, resulting that the two 
stars have si milar magnet osphere s, both wit h Alfven radius 
about 15 R. (|Trigilio et al . 2004; iLeto et al.: 2006). This re- 
inforces our decision to search evidence of CUVir-type co- 
herent emission from cr Ori E. 

The long history of the photometric, spectroscopic 
and magnetic variability presented by a Ori E is docu- 
mented i n the Catalog of observe d periods for the AP 
stars by Catalano fc R,ensonl lll984h and its supplements 
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[1991, 1993). Town send et al.l (|2010l ) have concluded that 
a Ori E is spinning down because of the magnetic braking 
and supplied the ephemeris: 

HJD = 2 442 778.829 4- 1.1908229£' + 1.44 x lO"^^^ [days] 

referred to the deeper light minimum. The spin down of 
a Ori E was already theor etically predicted by magnetohy- 
drodynamics simulations (|Ud-Doula. Owockv fc Townsendl 
I2OO9I ). 

The global magnetic fie ld of a Ori E is steady over three 
decades ( Oksala et al.l I2OIII ) and since iLandstreet fc Borral 
(| 19781 ) a pure magnetic dipole is assumed. Fig. [T] 
shows the magneti c field measurements of a Ori E by 
ILandstreet fc Borral (119781 ). iBohlender et all l|l987h and 
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Table 1. VLA observing log. 
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Figure 1. Fit of the literature magnetic field measurements - 
open circles. Legend of line style: single wave fit - dot line, double 
wave fit — dashed line. The solid horizontal line coincides with the 
null effective surface magnetic field. 



lOksala et al.l (|2011f ) phased with the previous ephemeris. A 
data fit with a single wave equation: 

B^ = Bo+ Bi sin 27r($ - $o) 

gives a reduced x^ = 12-4. Here Bo = 0.55 ± 0.02 [kG], 
Bi = 2.24 ± 0.03 [kG] and $o = 0.478 ± 0.002. 

From long time it is known that the variability of MCP 
stars can be accurately modeled through a sinusoidal wave 
and its first harmonic (jCatalano. Kroll fc Leonelll99in . In- 
deed a fit of magnetic data with the equation: 

Be = Bo + Bi sin27r($ - $o) + B2 sin27r(2$ - $1) 

results in a reduced x^ — 1.6, suggesting the non-dipolar 
nature of the aOriE magnetic field. Here Bo = 0.62 ± 0.02 
[kG], Bi = 2.21 ± 0.03 [kG], B2 = 0.63 ± 0.03 [kG], $0 = 
0.47 3 ± 0.002 and $1 = 0.563 ± 0.007. 

iLandolfi. Bagnulo fc Landi Degl'Innocentil (|l993) 

shown that the Bo variation due to a pure quadrupole 
is about 10% of the variability due to a dipole of equal 
polar strength. This means that a Ori E would present a 
quadrupole component at least two times stronger than the 
dipole one. 

The presence of multipolar magnetic components in 
(J Ori E is testified by the asymmetric emission features 
on the Ha wings. These periodically variable features 
have been ascribed to two circumstellar plasma clouds 
trapped in the magnetosph ere and co-rotating with the star 
(jLandstreet fc Borralll978l ): tra pped material was also rec- 
ognized from UV observations (jSmith fc Grootell2001r ). To 
explai n why one of the two f eature is stronger than the 
other, iGroote fc Hungen ([1983) suggested an asymmetry in 
the two clouds. In the case of 5 0riC, an MCP star twin 
of a Ori E, also presenti ng similar asymme tric emission fea- 
tures on the Ha wings, iLeone et al.l (|2010l ) mapped the cir- 
cumstellar matter distribution and ascribed the asymmetry 
between the two clouds to the magnetic quadrupole compo- 
nent necessary to explain the Stokes V profiles of the spectral 
lines. 

We conclude that the a Ori E phenomenology cannot 
be modeled as a simple dipole and that higher order com- 
ponents are necessary to explain the observational data. 



Code Frequencies 

[GHz] 



Epoeh 



AL346 5/15/43 

AT233 1.4/5 

AL56S 8.4/15 



95-Apr D 

99-Oct AB 

02-May AB 



3C48 0541-056 

3C48 0541-056 

3C286 0541-056 



AL618 


5/15/22/43 


04-Jan 


BC 


3C286 


0541-056 


Archive 


VLA data 










AL267 


5/8.4/15 


92-Oct 


A 


3C286 


0541-056 


AL348 


5/8.4/15/22/43 


95-Mar 


D 


3C286 


05324-075 


AL372 


5/8.4/15 


96-Mar 


C 


3C48 


0541-056 



3 RADIO OBSERVATIONS AND DATA 
REDUCTION 

Multi-frequency observations of a Ori E were carried out 
with the VLAJj in different epochs. Table [l] reports the in- 
strumental and observational details for each data set. 

To avoid significative phase fluctuations, scans on 
a Ori E were shorter than the Earth atmosphere coherence 
time and embedded between phase calibrator measurements. 
The 22 and 43 GHz observations were carried out using the 
fast switching mode between source and phase calibrator. 

The data were calibrated and mapped using the stan- 
dard procedures of the Astronomical Image Processing Sys- 
tem (AIPS). The flux density for the Stokes I parameter was 
obtained by fitting a two-dimensional gaussian (JMFIT) at 
the source position in the cleaned maps integrated over con- 
tiguous scans, the integration time ranges from about 10 
minutes to about 1 hour. As the uncertainty in the flux den- 
sity measurements we assume the r.m.s. of the map. 

The fraction of the circularly polarized flux density 
(Stokes V parameter) was instead determined performing 
the direct Fourier transform of the visibilities at the source 
position (DFTPL) without any temporal average. We are 
justified in this procedure by the absence of other circularly 
polarized sources in the field. Stokes V were later averaged 
with the same integration time of Stokes I. All VLA mea- 
surements are listed in Tabled] 



4 RADIO PROPERTIES OF a ORI E 

4.1 Light curves 

Fig. [2] shows the radio flux versus the rotational phase for 
!^ ^ 15 GHz. For completeness, we plot also the l itera- 
ture measuremen t s at 5 GHz by iDrake et al.l l|l987f ) and 
iLeone fc Umanal (119931). The top panel shows the vari- 
ability of Bp JOksala et al.ll201ll : iLandstreet fc Borra|[l978l : 

iBohlender et al.lll987l ). " 

As already kno wn (|Leonelll99ll : iLeone fc Umanalll993l : 
iTrigilio et akl l2004h the fight curve at 5 GHz of a Ori E 
is characterized by the presence of two maxima associated 
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Table 2. New and archival VLA observations. 



JD 

2400000+ 



Sl(Sv) 
m.Ty 



.JD 

2400000+ 



Sl(Sv) 
mjy 



njy 



JD 

2400000+ 



Sl(Sv) 
mJy 



JD 

2400000+ 



Sl(Sv) 
m.Jy 



nJy 



JD Si(Sv) 

2400000+ mJy 



nJy 



1^ = 1.4 GHz 



51453.833° 


2.60(+0.10 


0.15 


51453.862° 


2.50(+0.10 


0.15 


51453.896° 


2.10(+0.00 


0.15 


51453.951° 


2.00(+0.00 


0.20 


51453.998° 


1.70(+0.10 


0.20 


51454.050° 


1.80(-0.10 


0.20 


51454.135° 


2.30(+0.00 


0.20 


51455.827° 


2.00(-0.10 


0.15 


51455.862° 


2.30(+0.10 


0.10 


51455.902° 


2.40(+0.00 


0.10 


51455.941° 


2.80(-0.10 


0.10 


51455.998° 


2.70(+0.10 


0.10 


51456.067° 


2.70(+0.10 


0.10 


51456.120° 


2.70(+0.00 


0.10 


51456.162° 


2.80(+0.00 


0.10 


51458.841° 


2.00(-0.10 


0.08 


51458.899° 


2.20(+0.00 


0.10 


51458.956° 


2.30(-0.06 


0.07 


51459.031° 


2.60(+0.04 


0.08 


51459.122° 


2.40(+0.02 


0.08 


51461.929° 


2.70(+0.00 


0.10 


51461.976° 


2.80(+0.20 


0.10 


51462.021° 


2.80(+0.10 


0.10 


51462.067° 


2.80(+0.00 


0.10 


V 


=5 GHz 




48906.821" 


3.20(+0.00 


0.10 


48906. 973» 


4.10(+0.20 


0.10 


48907.125* 


4.90(-0.10 


0.10 


48908. 816» 


4.70(+0.20 


0.10 


48908.968'' 


3.20(+0.00 


0.10 


48909. 119°^ 


3.40(+0.10 


0.10 


48910. 811°^ 


3.30(+0.00 


0.15 


48910.962* 


3.60(+0.10 


0.15 


48911.114°^ 


4.30(+0.20 


0.15 


49803.498'' 


3.00(-0.20 


0.20 


49803.647*' 


3.80(-0.40 


0.20 


49827.591° 


4.50(+0.10 


0.10 


49833.450° 


4.30(-0.10 


0.10 


49837.398° 


3.00(+0.00 


0.10 


50170.328'' 


3.30(+0.00 


0.10 


50170.387'' 


3.97(-0.04 


0.08 


50170.446'' 


4.05(-0.10 


0.08 


50170.506'' 


3.93(-0.02 


0.07 


50170.565'^ 


4.28(-0.08 


0.07 


50170.624'' 


4.19(-0.02 


0.08 


50170.684'' 


3.78(-0.06 


0.08 


50171.380'' 


3.14(+0.03 


0.08 


50171.439'' 


2.90(+0.02 


0.08 


50171.499'' 


3.14(-0.02 


0.07 


50171.558'' 


3.55(+0.03 


0.08 


50171.617'' 


4.00(-0.06 


0.08 



50171.677'' 


4.10 


50172.377'' 


4.10 


50172.437'' 


3.43 


50172.496'^ 


3.14 


50172.555'' 


3.44 


50172.615'' 


3.31 


50172.674'' 


3.44 


51453.831° 


3.10 


51453.855° 


3.20 


51453.878° 


3.10 


51453.902° 


3.40 


51453.934° 


3.11 


51453.957° 


3.09 


51453.981° 


3.16 


51454.004° 


3.12 


51454.032° 


3.25 


51454.060° 


3.28 


51454.083° 


3.22 


51454.107° 


3.44 


51454.130° 


3.78 


51454.154° 


3.96 


51454.176° 


3.87 


51455.821° 


3.30 


51455.845° 


3.41 


51455.868° 


3.75 


51455.892° 


3.80 


51455.923° 


4.03 


51455.947° 


4.13 


51455.971° 


4.04 


51455.994° 


4.09 


51456.022° 


3.71 


51456.050° 


3.49 


51456.073° 


3.42 


51456.097° 


3.81 


51456.120° 


3.32 


51456.144° 


3.55 


51458.818° 


2.84 


51458.841° 


3.24 


51458.865° 


3.59 


51458.889° 


3.73 


51458.920° 


3.96 


51458.944° 


4.05 


51458.967° 


4.11 


51458.991° 


4.46 


51459.019° 


4.59 


51459.046° 


4.95 


51459.070° 


4.98 


51459.094° 


4.60 


51459.117° 


4.97 


51459.141° 


4.52 


51459.163° 


3.60 


51461.935° 


4.18 



(-0 


03) 


( + 


30) 


( + 


22) 


( + 


05) 


( + 


04) 


(-0 


01) 


( + 


OS) 


( + 


10) 


( + 


20) 
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10) 


( + 
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(-0 


02) 


( + 
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13) 


( + 


01) 


(-0 


07) 


(-0 


05) 


( + 


01) 


(-0 


05) 


( + 


00) 


(-0 


04) 


(-0 


13) 


( + 


20) 


( + 


15) 


( + 


01) 


( + 


20) 


( + 


IS) 


( + 


14) 


( + 


IS) 


( + 


15) 


( + 


17) 


( + 


26) 


( + 


33) 


( + 


20) 


( + 


06) 


( + 


17) 


(-0 


01) 


( + 


06) 


( + 


06) 


(-0 


07) 


(-0 


10) 


(-0 


06) 


(-0 


16) 


(-0 


11) 


(-0 


10) 


( + 


03) 


(-0 


06) 


(-0 


15) 


(-0 


OS) 


( + 


04) 


( + 


00) 


( + 


20) 



0.07 
0.10 
0.08 
0.08 
0.08 
0.07 
0.07 
0.10 
0.20 
0.10 
0.10 
0.08 
0.07 
0.07 
0.06 
0.05 
0.06 
0.06 
0.06 
0.06 
0.07 
0.08 
0.10 
0.08 
0.08 
0.10 
0.08 
0.08 
0.07 
0.08 
0.08 
0.07 
0.07 
0.08 
0.08 
0.08 
0.08 
0.07 
0.07 
0.07 
0.08 
0.07 
0.07 
0.07 
0.06 
0.06 
0.06 
0.07 
0.07 
0.07 
0.10 
0.06 



51461.958° 


4.42(+0.13 


0.05 


51461.982° 


4.14(+0.20 


0.06 


51462.009° 


4.14(+0.15 


0.06 


51462.037° 


4.21(+0.35 


0.06 


51462.061° 


4.00(+0.24 


0.05 


51462.085° 


3.88(+0.15 


0.06 


51462.111° 


3.53(+0.27 


0.05 


53031. 789^ 


3.40(+0.05 


O.OS 


u— 


8.4 GHz 




48906.869" 


4.40(+0.10 


0.10 


48906.887' 


4.40(+0.00 


0.10 


48906. 901» 


4.20(-0.30 


0.10 


48906. 915» 


4.00(-0.10 


0.10 


48906.934" 


3.90(-0.20 


0.10 


48907.011" 


4.30(-0.20 


0.10 


48907.030" 


4.60(-0.10 


0.10 


48907.048" 


4.70(^0.10 


0.10 


48907.067" 


4.80(+0.00 


0.10 


48907.086" 


4.90(+0.20 


0.10 


48907.163" 


5.00(+0.10 


0.10 


48908.854" 


3.90(+0.20 


0.10 


48908.873" 


3.70(+0.10 


0.10 


48908.891" 


3.50(+0.40 


0.10 


48908.910" 


3.30(+0.20 


0.10 


48908.929" 


3.20(+0.40 


0.10 


48909.006" 
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with the extrema of the magnetic field. We observe a similar 
modulation in the flux curve at 1.4 GHz. At both 1.4 and 
5 GHz, the minima in the light curves are close in phase 
with the null magnetic field, that occurs when the axis of 
the dipole is perpendicular to the line of sight. 

The shape of the light curves at 8.4 and 15 GHz are 



more complex and any relation with the Be variability is 
no longer simple. At 8.4 GHz more than two maxima are 
detected, two of these are clearly associated with the two 
maxima observed at $ « 0.3 and 0.95 in the 15 GHz light- 
curve. The amplitude of the light curves increases with the 
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Figure 2. VLA measurements of cr OriE plotted versus tiie ro- 
tation pliase *. Legend of data: AL267 (D), AL348 (v), AL346 
(t), AL372 (O), AT233 (■), AL568 (♦) and AL618 (a). Liter- 
ature data at 5 GHz are indicated as (-I-). Symbol size is equal 
to the flux error. Top panel reports the literature magnetic field 



frequency. It ranges from about 49% (at 1.4 GHz) to about 
82% (at 15 GHz) of the median. 

With the aim to detect CU Vir type pulses in a Ori E, we 
focus our attention in the phase windows where we expect to 
see them. Those windows are around the rotational phases 
corresponding to the null of the magnetic field curve. The 
first null is at $ Ri 0.9 — 0.1, the second at <l> « 0.3 — 0.5, both 
with a duration A<l>nuii ~ 0.2. The phase coverage of the ra- 
dio measurements around the first null is almost continuous 
in the frequency range 1.4 - 15 GHz, while around the sec- 
ond null there are fewer data. Assuming a pulse extends 
over a phase window A$ecm ~ 0.05, as for the narrowest 
pulse observed in CUVir, we estimate that the probability 
to have missed an ECM pulse in the first null is zero, since 
the maximum phase gap A<l>gap between contiguous points is 
smaller than A<1?ecm for all the observed frequencies. On the 
contrary, around the second null there are several gaps for 
which A$gap > A$ECM- In this case, we estimate the prob- 
ability to miss the ECM pulse as (A$gap — A$ECM)/A'I>nuii. 
Considering all the gaps in this phase window, those proba- 
bilities are 63%, 25%, 16% and 14%, respectively at 1.5, 5, 



8.4 and 15 GHz. For the whole rotational period, instead, 
we get 13%, 5%, 12% and 15%. 



4.2 Radio spectrum 

The 22 and 43 GHz VLA measurements of a Ori E do not 
cover the whole rotational period. However, our 22 and 43 
GHz observations are in coincidence with the steep flux in- 
crement observed at 15 GHz ($ « 0.95). We can then trace 
the complete radio spectrum from 1.4 to 43 GHz in the range 
of phases: 0.9 - 1, when the flux at 15 GHz approximately 
doubles. To minimize the effect of the lack of simultaneity of 
higher frequency measurements, we have fixed the phases: 
$ =0.867, 0.910, 0.922, 0.934 and 0.998, coinciding with the 
measurements at 43 GHz and with the interpolated value be- 
tween the two closest in phases 43 GHz observations. Fluxes 
at the other frequencies have been obtained by fits of the 
light curves, see Fig. (3] left panels. The five spectra are la- 
belled: (a), (b), (c), (d), and (e) in Fig. g] 

The radio spectrum (a), obtained before the rising 
phase of the 15 GHz flux, is peaked at intermediate frequen- 
cies (5 - 8.4 GHz). The radio spectra ((b), (c) and (d)) just 
before the 15 GHz maximum are, instead, peaked at higher 
frequency (15 GHz) and show a steep rise at 43 GHz (Fig.|3l 
right top panel). At frequencies below 15 GHz the spectral 
index is positive, with a ~ 0.3 (5*^ oc v"'). The spectrum (e), 
obtained after the maximum, closely resembles the spectrum 
(d) (having similar 15 GHz flux level but obtained before the 
maximum) at the lowest frequencies, whereas at the highest 
frequencies (22 and 43 GHz) fluxes are decreasing. 



4.3 Circular polarization 

The circularly polarized radio emission shows two extrema of 
opposite signs (Fig. [l] bottom panels), detected above the 
3(7 detection threshold, associated with the magnetic field 
extrema (Fig. [1} . The degree of polarization enhances as the 
frequency increases. Positive degree of circular polarization 
is detected when the north magnetic pole is close to the line 
of sight and negative in presence of the south pole. When 
the magnetic poles are close to the direction of the line of 
sight we observe most of the radially oriented field lines. In 
this case the gyrosynchrotron mechanism give rise to radio 
emission partially polarized, respectively right hand for the 
north pole and left hand for the south pole. 

The overall behavior of the circularly polarized radia- 
tion from a Ori E r esembles so close ly the case of CU Vir at 
5, 8.4 and 15 GHz (|Leto et aLlbOOeP . This can be considered 
the typical behavior of the gyrosynchrotron emission from a 
magnetosphere characterized by a mainly dipolar symmetry. 



5 EFFECT OF THE MAGNETIC 

QUADRUPOLE ON THE RADIO EMISSION 

5.1 Comparison with the 3D model 

In order to investigate the radio emission from MCP stars, 
in lTrigilio et al.l l|2004 ) we developed a 3D model to compute 
the gyrosynchrotron emission from a magnetosphere shaped 
in the framework of the o blique rotator mo del, that assumes 
a simple magnetic dipole. lLeto et al.l (|2006r ) extended the 3D 
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Figure 3. Left panels: zoom of the radio light curves at phases where high frequencies data are available; vertical lines (a, b, c, d, e) 
show the phases where the spectra have been realized. Right panels: radio spectra at phases (a), (b), (c), (d), and (e); rising phase up 
panel, decay phase down panel. 



model to solve the radiation transfer for the circularly po- 
larized emission (Stokes V parameter). Such a 3D modeling 
has been suc cessfully applied to a Ori E and H P 37017 by 
iTrigilio et al] (i2004 1 and to CU Vir bv lLeto et al. (2006 ). 

In this model a radiatively driven stellar wind flows 
along the magnetic field lines. The wind channeled by the 
magnetic field reaches the magnetic equator. At the Alfven 
radius, where the kinetic energy density of the wind exceeds 
the magnetic one, the stellar wind breaks the magnetic field 
lines generating current sheets, where particle acceleration 
occurs. The energetic electrons moving towards the stellar 
surface emit by gyrosynchrotron mechanism. This thin layer 
is named middle magnetosphere and it is connected to the 
annular regions around the magnetic poles at high magnetic 
latitudes. At lower magnetic latitudes the wind is instead 
confined in a dead zone (inner magnetosphere) , resulting in 
an accumulation of circumstellar matter whose density de- 
creases outward, i n accord with the rnagnet ically confined 
wind shock model (jBabel fc Montmerlelll997^ . This circum- 
stellar material absorbs the continuous gyrosynchrotron ra- 
diation, resulting in the observed deep rotational modula- 
tion. 

Here we adop t the model parameters given by 



ITrigilio et akl (|2004l ). used to fit the literature 5 GHz Ught 
curve of a Ori E, to compute the expected radio fiux den- 
sity and fraction of circular polarization (ttc = Sv/Si) for 
all the frequencies. In Fig. [4] we have compared the simula- 



tions with the observations. The 1.4 GHz Stokes I variability 
is reproduced in shape, however slightly overestimating the 
average fiux. As to the highest frequencies, our 3D model 
fails to reproduce the observed Stokes I variability. On the 
contrary, Stokes V simulations are in good agreement with 
the measurements. 

Since the source is optically thick, the radiation at dif- 
ferent frequencies probes the stellar magnetosphere at differ- 
ent depths, where r « 1. Higher frequencies are mainly emit- 
ted in the inner regions, where magnetic field is stronger. 
In particu lar, inspecting the m odeled radio maps of a Ori E 
shown bv ITrigilio et all l|2004l) . the bulk of the 5 GHz gy- 
rosynchrotron emission arise from magnetospheric regions 3 
- 8 R* away from the stellar surface, whereas the 15 GHz 
originates from inner layers (1 - 5 R«). 

We conclude that the 3D model reproduces the behavior 
of the 1.4 and 5 GHz radio emission, coming mainly from 
regions far from the star, where the dipolar component of 
the magnetic field dominates. The nearly-dipolar symmetry 
can be also identified from the behavior of the polarized 
emission, well reproduced by the model simulations at all 
the frequencies (Fig. |3| bottom panels). At 8.4 and 15 GHz 
the model instead fails to reproduce the observed Stokes 
I modulation, as expected if a non dipolar magnetic field 
component is present. 

In the case of gyrosynchrotron mechanism, the radio 
emission from a region with almost parallel magnetic field 
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Fi gure 4. Comparison of the observed radio light-curves of cr OriE and the theoretical curves (dashed line) computed with the 3D model 
bv TTrigilio et alj l|2004l'l . 



lines directed toward us (North), or in opposite direction 
(South), is right-hand (RCP), or left-hand (LCP), circularly 
polarized; Stokes V, which is the difference between RCP 
and LCP, is positive, or negative; Stokes I, which is the sum 
of RCP and LCP is sensitive to both polarizations. In a 
region with dishomogeneous magnetic field on a small scale, 
the polarization is, in average, null. When a quadrupolar 
field is present, we see always regions with antiparallel field 
lines; Stokes V is always null, while Stokes I depends on 
the strength of each individual component. In this simple 
way we can explain the behavior of Stokes V, which well 
agrees with the dipolar model at any frequency, and of the 
low frequency Stokes I, arising from regions where only the 
dipolar field is relevant. The high frequency Stokes I, instead, 
is sensitive to the quadrupolar component, important in the 
low magnetosphere. We are supported in such a conclusion 
from the Stokes I and V variability p resented at 5 GHz by 
HR5624 (jLim. Drake fc Linskvlll996|). that is cha racterized 
by a multi-polar magnetic field ( Landstreetlll99d ). 

The spectra shown in Fig. [3] evidence as the radio flux 
density at the higher frequencies steeply rise going up to the 
maximum at phase ~ 0.95, clearly observed in the 15 GHz 
light curve. This could be a consequence of strong dishomo- 
geneity in the magnetic field close to the star. 

5.2 The absence of the cyclotron maser 

This work is aimed mainly at searching for CUVir-like 
pulses from crOriE. In the following we will compute the 
expected frequency range of the ECM in our target. 

The minimum maser frequenc y observed in CU Vir is 
610 MHz jStevens fc Georgel [2OI0I '). the maximum is 2.5 
GHz. Since the maser frequency is a harmonic of the lo- 
cal gyrofrequency, v = s x vb, with s = 2, and assuming 
a simple dipolar fleld, for which the magnetic field strength 
above the pole is given by 

Bdip = Bp(R*/-R)^ 

we can deduce that for CUVir (_Bp = 3000 G) the maser 
at 610 MHz is generated at 2 R* above the stellar surface. 
The shaded area of Fig. [S] defines the frequencies and the 
corresponding distances associated with the ECM pulses as 
observed on CUVir. At the same height, for aOriE, with 



a polar field strength Bp = 6800 G, the corresponding fre- 
quency of emission is 1.4 GHz. Therefore the cyclotron maser 
emission from o OriE, if it does exist, should be observable 
at frequencies higher than 1.4 GHz. The maximum possible 
frequency is given by the second harmonic of the gyrofre- 
quency very close to the stellar surface, corresponding to 
about 38 GHz. 

Moreover, a necessary condition for the ECM is that 
the local plasma frequenc y must be lower than t he local gy- 
rofrequency, i.e. Vp <C vb IjMelrose fc Dulklll982l ). Since the 
plasma frequency is given by t-p = 9 x 10~^ne GHz, the 
maser emission near to the magnetic pole is possible only 
if the local thermal electron number density is less than 
about 4 X lO^'^ cm""'. In the current model the thin cav- 
ity where the maser should originate is the middle magne- 
tosphere above the magnetic pole, where the local plasma 
is the ionized therm al wind that flows a long the magnetic 
field lines. Following iTrigilio et al.l l|2008r ) and adopting the 
characterist ic values of the wind derived with the 3D model 
for cr OriE (|Trigilio et al.l 120041 ) ■ we estimate that the den- 
sity of the wind is about lO'^^ cm~'^ near the stellar surface, 
well below the theoretical quenching limit. As the magnetic 
field strength decreases going outward, the plasma frequency 
equates the gyrofrequency at approximately 7R, above the 
surface. Therefore, we can conclude that the condition for 
generation of the cyclotron maser is satisfied in the middle 
magnetosphere of o Ori E. 

The radial behavior of the plasma density in the mid- 
dle and inner magn etosphere, is s hown in Fig. [5] for a Ori E 
dTrigilioet a"Lll2004l ) and CU Vir (iLeto et al.ll2006[). As re- 
centl y observed on CUVir l|TrigiUo et al.l 1201 ll ; ILo et al] 
[2013) the ray path of the ECM, which is polarized in the 
x— mode, is refracted by the plasma trapped in the inner 
magnetosphere, the refractive index is given by rirofr ~ 
^JY^^l^J2v^, with i/'ECM ~ 2vb- Since the conditions of 
the plasma are quite similar for a Ori E and CU Vir (Fig. [5]), 
we expect an upwards refraction. In few words the radiation 
should be deviated, not absorbed. 

Despite the light curves in the expected frequency range 
of the ECM (between 1.4 and 15 GHz) are well sampled, in 
particular near the rotational phases where we expect to de- 
tect it, i.e. in coincidence with at least one null of the mag- 
netic curve (Fig. [2]), we do not see any hint of the pulses as 
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Figure 5. Radial dependence of the thermal electron density 
number inside the "middle" and the "inner" magnetosphere, in 
the cases of crOriE and CU Vir. In the top x-axis the ECM fre- 
quencies {2v^) for the two stars, are shown, the shaded area iden- 
tifies the corresponding magnetospheric regions. 



observed on CU Vir. Since the frequency of 1.4 and 15 GHz 
corresponds to the second harmonic of the local gyrofre- 
quency at 2 and 0.4 R« above the pole respectively, we can 
conclude that in this region the conditions needed to trigger 
the cyclotron maser are not fulfilled. 

As highlighted in this paper, the magnetic field of 
(7 Ori E is characterized by the presence of a quadrupolar 
component, whose strength Squad, which is prevailing in 
the stellar surface, decreases faster than the dipolar com- 
ponent (Bquad oc i?~^, Bdip oc R~^). Therefore the actual 
magnetic vector in the low magnetosphere can be quite dif- 
ferent from the simple dipolar one. As already discussed, this 
quadrupolar component significantly affects the gyrosyn- 
chrotron emission at high frequency, originating in the deep 
magnetospheric layers. In a similar way, the quadrupolar 
component of the magnetic field, stronger then the dipolar 
one at low altitude, affect the propagation of the electrons 
moving towards the stellar surface, causing the moving up 
of the magnetic mirroring point and consequently different 
interaction with the thermal plasma of the stellar magne- 
tosphere, and possible changing of the ECM conditions. In 
fact, for a charged particle moving with a pitch angle in a 
not uniform magnetic field, its magnetic moment is invari- 
ant, i.e. sin"^ <j)/B — constant, leading to magnetic mirroring 
when the particle moves in a converging magnetic field. As- 
suming for G Ori E a simple dipolar field, a non-thermal elec- 
tron injected with <jf)o ~ 1° at the Alfven radius (i?Aitvcn ~ 15 
Rt), where BAifven ~ 2 G, will reach the stellar surface with 
a pitch angle (f> ~ 90° (sin^ (f> = sin^ 4>o'xBp/BAiivcn)- If there 
is a magnetic quadrupolar component, which is two or three 
times stronger than the strength of the dipolar field, the 
same electron is refiected at 0.3 - 0.35 R* from the surface, 
and cannot reach the deeper and denser magnetospheric lay- 
ers; in this case it is refiected back without being absorbed. 
Only those electrons injected at i?Aifven with a very small 
pitch angle (^ 1°) can reach the stellar surface. Accord- 
ingly, the fraction of non-thermal electrons that can develop 
an anisotropic pitch angle distribution is very small, and 



then the loss cone angle is already closed when the electrons 
reach the region where the maser at frequency below 15 GHz 
should arise. The quadrupole component vanishes 5 - 6 R, 
away from the surface and this explain the observed a Ori E 
behavior. 



6 CAN O-TYPE STARS EMIT ECM PULSES? 

Hot 0-type stars have strong ionized stellar winds and , 
in few cases, organized magnetic fields (|Donati et al.ll2002l . 
l2006bl ). like MCP stars. Such characteristics imply that the 
magnetic 0-type stars may be sources of ECM. 

Typically MCP stars have weak winds (M « 10"^° - 
10~^ M© yr~^) strongly dominated by the magnetic fields 
{Bp ~ 1-10 kG), with large co- rotating magnetospheres. 
The measured magnetic field in t he 0-type stars is about 
1 kG or less ( Hubrig c t al.l l201ll ) and the mass loss rate 
is higher then 10"^ M p yr~\ Since i^Aifvcn oc B^'^/M'^^^ 
(|Ud-Doula et al.ll2009l '). the very massive wind and the rel- 
atively weak magnetic field locates the Alfven surface of 
the known magnetic O-type stars close to the stellar sur- 
face (< 2 R*). In the case of single magnetic 0-type 
stars, shocks originating in the wind channeled by the mag- 
netic field can accelerate electrons up to relativistic energy 
(|Van Loo. Runacres fc B lominc 200j|). These, in turn, could 
emit at the radio regime by gyrosynchrotron mechanism. 
On the other hand the 0-type stars are characterized by 
radio emission at centimeter wavelengths mainly ascribed 
to the thermal free-free emission from the ionized stellar 
wind, whose typical spectral index is a « 0.6 (Si, oc i/") 
(I Wright fc Barlow! Il975h . The radio photosphere, at cen- 
timet er wavelengths, has a radius of « 100 R* (JBlommg 
l201l[ ). meaning that the plasma of the wind is optically 
thick inside and that any non-therm al emissio n generated 
within is completely absorbed (Van Loo et al.l |2006;) . The 
ECM from O-type stars, if any, would suffer from the same 
absorption effect. 

Moreover, for a main sequence magnetic O-type star 
with stellar radius Rt — 7 — 8 Kq, radiative wind mass 
loss rate M = 10"® M© yr"\ terminal velocity Uoo = 2000 
km s~^ and dipolar magnetic field with Bp = 1000 G, the 
condition Vp > vs, is valid everywhere. This is a further con- 
sideration against the possibility that a magnetic O-type 
star can emit ECM pulse. Up to now the signatures of non- 
thermal emission mechanism, flat or negative spectral index 
and variability, in the O-type stars are observed only in mul- 
tiple systems, the relativistic electrons origin is thus ascribed 
to the acceleration in the shocks that ta ke place in the col- 
liding winds far from the stellar surface l|Blommell2011f ). 

In conclusion, the conditions for generation and propa- 
gation of the ECM emission are not satisfied for the known 
magnetic O-type stars. The behavior of the magnetized O- 
type stars is thus significantly different from the MCP stars 
(B/A spectral type) in the radio domain. The key to ex- 
plain when the transition occurs is probably the shaping of 
the stellar magnetosphere, which depends on mass loss rate, 
magnetic field strength and rotational period. 
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T CONCLUSIONS AND OUTLOOK 

CUVir is up to now the only known magnetic chemically 
peculiar star presenting electron cyclotron maser emission. 
The characteristics of the MCP star a Ori E suggest that it 
should present the same phenomenology of CUVir at fre- 
quencies higher than 1.4 GHz. We have obtained multi fre- 
quencies radio observations during the whole the rotation 
period of a Ori E, in particular at the phases where the ECM 
is expected. However, no cyclotron maser emission has been 
detected in this star. 

There are indications that the magnetic field topology 
of (J Ori E cannot be considered as a simple dipole. The pres- 
ence of the quadrupole field could inhibit the development 
of the conditions able to power the cyclotron maser emis- 
sion. Nevertheless, we cannot rule out the possibility that 
the ECM could occur at other frequencies not analyzed in 
this paper, such as the 2.5 GHz, or at phases near the second 
null of the magnetic curve that are not fully covered. 

On the other hand, to better understand the ECM 
in the wider context of plasma processes it will be ex- 
tremely interesting to extend such investigation to a big- 
ger sample of MCPs. The average radio lumi nosity of the 
MCP stars is about 10^'^-^=^°-^ Jergs s~^ Hz'^] (JDrake et all 
ll987l : lLinskv. Drake fc Bastianlll992l ): it has been also shown 
that the radio luminosity o f the MCP stars increases with 
the effective tem perature (jLinskv. Drake fc BastianI 11993 : 
iLeone et al.l 1 1994) . The radio interferometers of new gen- 
eration, like the EVLA or the forthcoming Australian SKA 
Pathfinder (ASKAP, which wiU operate at 1.4 GHz), will 
allow to reach the detection limit of few /ijy in all sky deep 
surveys. Assuming a threshold of 10 ^Jy we estimate in 
about 2000 pc the maximum distance within will be pos- 
sible to detect the radio em ission from the MCPs. Following 
iRenson fc ManfroidI ((20091) we can assume the MCP stars 
uniformly distributed in space. It is reasonable to expect 
that the number of radio detections will increase of about 
an order of magnitude, giving the opportunity to get a larger 
statistics of the physical conditions of the magnetospheres, 
like the magnetic field strength and orientation and the ther- 
mal plasma density, to correlate with the ECM. 

We want to stress here that this class of objects provides 
the unique possibility to study plasma process in stable mag- 
netic structures, whose topologies are quite o ften well deter- 
mine d by several independent diagnostics IjBvchkov et al.l 
I2OO5I ). thus overcoming the variability of the magnetic field 
that is one of the major problem that prevents an accu- 
rate modeling of ECM in very active stars such dMe or 
close binary systems. On the other side, moving towards 
earlier spectral types, stronger and denser stellar winds to- 
gether with a weaker magnetic field should inhibit the onset 
of the ECM instability. Moreover, as shown for the proto- 
type CUVir, the observations of persistent coherent pulses 
of ECM, which act as a clock of the star, in other MCP 
stars would offer a valuable tool for the angular momentum 
evolution of the objects belonging to this class. 
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